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Abstract 



We study the motion-planning problem for a car-like robot whose turning radius is bounded from 
below by one and which is allowed to move in the forward direction only (Dubins car). For two robot 
configurations a, a' , let l(a, a 1 ) be the shortest bounded-curvature path from a to a' . For d ^ 0, let 
£(d) be the supremum of £(a, a'), over all pairs (a, a') that are at Euclidean distance d. We study 
the function dub(d) = £(d) — d, which expresses the difference between the bounded-curvature path 
length and the Euclidean distance of its endpoints. We show that dub(d) decreases monotonically 
from dub(0) = 7n/3 to dub(d*) = 2tv, and is constant for d ^ d* . Here d* ~ 1.5874. We describe 
pairs of configurations that exhibit the worst-case of dub(d) for every distance d. 

1 Introduction 

Motion planning or path planning involves computing a feasible path, possibly optimal for some criterion 
such as time or length, of a robot moving among obstacles; see the book by Lavalle [17] and book chapters 
by Halperin et al. [13] and Sharir [23]. A robot generally comes with physical limitations, such as bounds 
on its velocity, acceleration or curvature. Such differential constraints restrict the geometry of the paths 
the robot can follow. In this setting, the goal of motion planning is to find a feasible (or optimal) path 
satisfying both global (obstacles) and local (differential) constraints if it exists. 

In this paper, we study the bounded- curvature motion planning problem which models a car-like robot. 
A car (with front-wheel steering) is constrained to move in the direction that the rear wheels are pointing, 
and it has a hxed maximum steering angle. This makes the car travel in a motion with fixed minimum 
turning radius, which means that the car must follow a curvature- constrained path. More precisely, we 
have the following robot model: 

Robot model (Dubins car). The robot is considered a rigid body that moves in the plane. A 
configuration of the robot is specified by both its location, a point in R 2 (typically, the midpoint of the 
rear axle), and its orientation, or direction of travel. The robot is constrained to move in the forward 
direction, and its turning radius is bounded from below by a positive constant, which can be assumed to 
be equal to one by scaling the space. In this context, the robot follows a bounded-curvature path, that is, 
a diffcrentiable curve whose curvature is constrained to be at most one almost everywhere. 

Planning the motion of a car-like robot has received considerable attention in the literature. In this 
paper, we consider the cost of this restriction: How much longer is the shortest path made by such a 
robot compared to the Euclidean distance travelled? 

Formally, consider two configurations a and a'. Let £(a,a') denote the length of a shortest curvature- 
constrained path from a to a', and let d(a,a') denote the Euclidean distance between a and a 1 . We 
define 



Note that the supremum here is not a maximum, as the path length is not a continuous function of the 
orientations at the two endpoints. Our goal is to understand the function dub : R i— >• K in detail. While 
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dub(e?) = sup{£(<7, a') — d | a, a configurations with G?(er, a') = d}. 



(1) 



this is a natural and fundamental question related to motion planning with bounded curvature, it is also 
a relevant question that has repeatedly appeared in the literature, with only partial answers so far. 

Dubins showed that the shortest curvature-constrained path between two configurations consists 
of at most three segments, each of which is either a straight segment or a circular arc of radius one. 
Using ideas from control theory, Boissonnat et al. [7], in parallel with Sussmann and Tang [35], gave 
an alternative proof. Sussmann [25] extended the characterization to the 3-dimensional case. Bui et 
al. [nil discussed how the types of optimal paths partition the configuration space, and also proved that 
optimal paths for free final orientation have at most two segments [6j. Significant work has been done on 
the problems of deciding whether a bounded-curvature path exists between given configurations among 
different kinds of obstacles and finding the shortest such path [I ^ I51 l2TllT51 l51l3|l21IT]|gilg]. 

At least two interesting problems have been studied where not configurations but only locations 
for the robot are given. The first problem considers a sequence of points in the plane, and asks for 
the shortest curvature-constrained path that visits the points in this sequence. In the second problem, 
the Dubins traveling salesman problem, the input is a set of points in the plane, and asks to find a 
shortest curvature-constrained path visiting all points. Both problems have been studied by researchers 
in the robotics community, giving heuristics and experimental results [2"2"l Hl?l |2"0"] . From a theoretical 
perspective, Lee et al. [THj gave a linear-time, constant-factor approximation algorithm for the first 
problem. No general approximation algorithms are known for the Dubins traveling salesman problem 
(the approximation factor of the known algorithms depends on the smallest distance between points) . 

All this work depends on some knowledge of the function dub. Lee et al. [IS], for instance, prove that 
the approximation ratio of their algorithms is max(_4, tt/2 + B/ir), where A = 1 + sup{dub(d)/d | d ^ 2} 
and B = sup{dub(d) + d | d 2}. They claim without proof that dub(d) ^ 2tt for d ^ 2 and derive from 
this that A = 1 + 7r and B ^ 5ir/2 + 3, leading to an approximation ratio of about 5.03. We give the first 
proof of A = 1 + 7r, and improve the second bound to B — 2 + 2tt, improving the approximation ratio of 
their algorithm to 2 + 2/ir + tt/2 « 4.21. 

Savla et al. [23] prove that dub(d) ^ kit, where k € [2.657, 2.658], and conjecture based on numerical 
experiments that the true bound is 77r/3. We show that this is indeed true. 

Results. We show that d dub(d) is a decreasing function with two breakpoints, at \/2 and at d* m 
1.5874 (see Figure [T|. More precisely, we have dub(0) = 77r/3 and the two breakpoint values are 
dub(v / 2) = 57r/2 — \72, and dub(d*) = 2tt. The function dub(d) is constant and equal to 2tx for d^? d* . 
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Figure 1: The graph of the function dub(d). 



For ^ d < \/2 and for d ^ d* , the supremum in ([T]) is in fact a maximum, and we give configura- 
tions (J, a' at distance d such that £(a, a 1 ) = dub(c?) + d. Perhaps surprisingly, for \[2 ^ d < d* , there are 
no such configurations — the supremum is not a maximum. 
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Figure 2: An example of a pair of configurations. 



2 Preliminaries 

Notations. For two points P and Q, we denote by PQ the line segment with endpoints P and Q, and 
by PQ an arc of unit radius with endpoints P and Q. (If the length of PQ is less than two then there 
are four such arcs, so unless it is clear from the context, we will specify the supporting circle and the 
orientation of the arc.) We denote the length of the segment PQ as \PQ\ or simply as \PQ\, and the 
length of the arc PQ as | PQ | . 

Without loss of generality, we assume that the starting configuration is (0, 0, a) — that is, we start 
at the origin S = (0, 0) with orientation a — and the final configuration is (d, 0, j3) — that is, we arrive at 
F = (d, 0) with orientation j3. Here, a and /3 express the orientation of the robot as an angle with the 
positive x-axis, and d ^ is the Euclidean distance of the two configurations. 

The open unit (radius) disks tangent to the starting and final configurations are denoted Ls, Rs, Lp, Rp, 
where the letters L or R depend on whether the disk is located on the left or right side of the direction 
vector (see Figure [2]). 

Let £sj r S^FT r F denote the centers of Ls, Rs, Lp, Rp, respectively. For future reference, we note 
their coordinates: 



is = (cos(a + 7r/2), sin(a + tt/2)) = (— sin a, cos a) (2) 

rs = (cos(a — 7r/2), sin(a — tt/2)) = (sin a, — cos a) (3) 

l F = (d + cos( / 3 + tt/2), sin(/3 + tt/2)) = (d - sin /3, cos/3) (4) 

r F = (d + cos(,3- tt/2), sin(/3- tt/2)) = (d + sin /?, - cos /?). (5) 



Distances between centers. The following distances will be frequently used: 
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It will be convenient to name the following terms: 

D L = d L ^l-(d L /4) 2 (10) 
D R = d nV /l - (rf R /4) 2 (11) 

Dubins paths. Dubins [TT] showed that for two given configurations in the plane, shortest bounded- 
curvature paths consist of arcs of unit radius circles (C-segments) and straight line segments (S'-segments); 
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Figure 3: (a) Monotonicity of the length function of CSC-paths; (b) Non-monotonicity of the length 
function of CCC-path. 

moreover, such shortest paths are of type CCC or CSC , or a substring thereof. These types of paths are 
referred to as Dubins paths. 

For given d ^ 0, a, j3 € [0, 2ir], there are up to six types of Dubins paths. The two path types LSL 
and RSR use outer tangents — these path types exist for any choice of d,a,/3. The two path types LSR 
and RSL use inner tangents, and exist only when the corresponding disks are disjoint. In particular, LSR 
exists if and only if d LR ^ 2, and RSL exists if and only if e? RL ^ 2. The remaining two path types LRL 
and RLR exist whenever there is a disk tangent to the two disks, and so LRL exists if and only if d L ^4, 
and RLR exists if and only if d R ^ 4. 

Dubins showed that in LRL- and i?Li?-paths the middle circular arc has length larger than ir. This 
implies that of the two unit radius disks tangent to L$ and Lp, only one is a candidate for the middle 
arc of an LRL-p&th, and similar for i?Li?-paths. 

For d > and ^ a, (3 ^ 2ir, we define lsl(c?, a, j3) to be the length of the LSL-path from (0, 0, a) 
to (d,0,/3). We define RSR, LSR, RSL, LRL, RLR similarly, defining the length to be oo if no path of that 
type exists. The length of the shortest bounded-curvature path from S to F is then 

£(d, a, (3) = min { LSL(d, a, fi), RSR(d, a, (3) , lsr(cZ, a, (3) , RSL(<i, a, j3), lrl(g£, a, (3) , rlr(c£, a, /?)}, 

and our goal is to bound dub(<i) = sup ^ Q p^in ^(^> a ) /^) — ^- (Note that the supremum here is not always 
a maximum as the function £ is not continuous.) 

We will often suppress the argument d for these functions when the distance d is fixed and understood. 

Monotonicity of the Dubins cost function. Let □ = [0, 2ir] 2 denote the range of (a, /3). Consider 
two distances d\ < di and (a, j3) G □, and assume that we have a bounded-curvature path from (0, 0, a) 
to (di, 0, /3) of length I ^ dub(di) +d\. If this path has a horizontal tangent where the path is oriented to 
the right (in the direction of the positive x-axis) , then we can insert a horizontal segment of length d 2 — d\ 
at this point, and obtain a path from (0, 0, a) to (e^i 0, /3) of length t + (rf 2 — di) ^ dub(di) + di. See, for 
instance, Figure ^ a). 

If this was possible for all (a, f3) E □, then it would imply that dub(c?2) ^ dub(di), and it would follow 
that the Dubins cost function is monotone. Unfortunately, not all Dubins paths have horizontal tangents 
with the correct orientation (see Figure |3^b) for an example), and so proving the monotonicity of the 
Dubins cost function will require much more work. However, we can start with the following lemma: 

Lemma 1. Let d\ < di, and {a, 13) € □. // there is a path of length £ of type RSR, LSL, LSR, or RSL 
from (0, 0, a) to (d\, 0, (3), then there is a path of length £ + {di — d{) from (0, 0, a) to (di, 0, /?). 

Proof. It suffices to show that any of these path types must have a horizontal tangent oriented in the 
positive x-direction. By symmetry, it suffices to show this for RSR- and RSL-p&ths. The topmost point 
on a .R.S'iJ-path necessarily has the correct tangent, so consider an RSL-paXh. It consists of a right-turning 
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(a) Original path (b) Mirrored along x-axis 




(c) Point symmetry in (tt, tt) 



Figure 4: Symmetry by mirroring the path along the x-axis. 

arc STi on i?s, a segment TiT 2 , and a left-turning arc T 2 F on Lp. If STi contains the topmost point 
of i?5, or if T2F contains the bottommost point of Lp, these have the correct tangent, and we are done. 
Otherwise the path cannot possibly reach the positive contradiction. □ 



Symmetries. For a fixed d > 0, determining dub(eQ essentially amounts to finding (a, (3) € □ maxi- 
mizing £(a,f3) ("essentially" since the maximum may not actually be assumed). We now observe that the 
function £(a, f3) has two symmetries. 

First, we can mirror a path around the x-axis. This maps a to —a, (3 to —j3, left disks to right disks, 
and right disks to left disks. As a result, we have, say, LSR(d, a, /?) = RSl(g?, ~a, — ft), and in general we 
have £(d, a, (3) = l(d, -a, -/?). See Figure [IJ 

Second, we can mirror a path around the line x = d/2 and reverse the direction of the path. If 
the original path connected (0, 0, a) with (d, 0, /3), the new path connects (0, 0, 2tt — (3) to (d, 0, 2ir — a). 
The transformation maps left disks to left disks and right disks to right disks, so we have, for instance 
LSR((i, a, f3) = LSR(d, 2tt — ft, 2ir — a), and in general £(d, a, (3) = £(d, 2it - p, 2ir — a). See Figure [5] 

Considered as symmetries on □, the mapping (a,/?) (—a, —ft) is a point symmetry in (tt, 7r), while 
the mapping (a, /3) t-¥ (2ir — f3, 2n — a) is a reflection around the line (3 = 2~k — a. 

It follows that sup( Qj ^) g [ 2Tr ] 2 ^(^, a, /3) = supr a p\ eA £(d, a, /3), where A is the triangle with corners 
(0,0), (tt, tt), and (0, 2tt), or in other words the region 

A : 0^a^7r and a < (3 sC 2ir - a. 



In the following we will thus be able to restrict our considerations to the triangle A (see Figure 6(a) ). 



A new parameterization. We now introduce a new parameterization of the (a, /3)-plane, which will 
sometimes be more convenient to work with: 

j3 + a . f3 — a 

c = — ^ — i = — ^ — • 

2 ' 2 

In other words, we have 

a = a — 6, j3 = a + S. 
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(a) Original path 



(b) Mirrored along the line x = d/2 while reversing 
the direction 




(2tt - /3, 2tt - a) 



a 2ir 

(c) Reflection around the line f} = 2n — a 

Figure 5: Symmetry by mirroring the path along the line x — d/2 while reversing the direction. 
Recall our triangle A from above. In the (cr, ^-representation, the triangle A is the triangle 

A : sC a < tt and ^ S ^ a, 



or the bottom right half of the square T = [0,7r] 2 (see Figure 6(b)). In this representation, our first 
symmetry maps (a, 6) to (— cr, — 5), while the second symmetry maps (a, 5) to (— cr, S). We thus have 
point symmetry in the origin, as well as mirror symmetry around the ^-axis. In addition, (cr + tt, 5 + tt) 
represents the same angles as (cr, 5) since a = a — S and f3 = a + 5, and so we also have point symmetry 
in the point (7r/2, 7r/2), or in other words £(d, cr, S) = i(d, it — a, tt — 5). 

Distances between centers (using a and 5). The following lemma will allow us to express the 
center distances in terms of cr and 5: 




J-axis 



a 

(a) (a, /^-representation 




cr TT 

(b) (<r, <5)-representation 



Figure 6: A in two different representations. 







Figure 7: For d = 0, a = and /3 = 7r, we have lrl(o!, a, /3) = RLR(d, a, ft) — ^p. 

Lemma 2. Let d J? 0, and define the points p($) = (cos $, sin - !?) and q{d) — (d — cos sin??) on the unit 
radius disks around (0,0) and (d, 0). Then 

|p(tf - + </>)| 2 = d 2 - 4dcos??cos</> + 4cos 2 tf. 

Proof. We have 

0)o(^ + 0)| 2 = (cos(0- 0) -d + cos(?9 + </>)) 2 + (sin(i?- 0) -sin^ + 0)) 2 
= (cos i? cos <f) + sin $ sin <fi — d + cos $ cos — sin i? sin </>) 2 

+ (sin ?? cos <j> — cos i? sin — sin cos ^ — cos $ sin </>) 2 
= (2 cos cos - d) 2 + A cos 2 ?9 sin 2 <j) 
= d — 4d cos i? cos + 4 cos $ cos <p + 4 cos i? sin 



= d 2 - 4o! cos i9 cos + 4 cos 2 d. □ 

Lemma [2] leads to the following expressions for the squared distances between our disk centers: 

dl = \e s ( F \ 2 = d 2 -4dsin<5coso- + 4sin 2 (5 (12) 

d 2 = |r s r F | 2 = d 2 + Ad sin S cos a + 4 sin 2 6 (13) 

din = \tsrF I 2 = d 2 + Ad cos S sin a + A cos 2 5 (14) 

d 2 L = \r s £F\ 2 — d 2 - AdcosSsina + 4cos 2 S (15) 

To see this, observe that is — p(a + vr/2), r$ = p(a — 7r/2), lp = q{n/2 — j3), and rp = q(3n/2 — (5). 



For d L , set = 7r/2 — (5, = — cr; for d R , set i9 = 7r/2 — S, = 7r — cr; for d LR , set $ = 7r — (5, <f) = tt/2 — <r; 
for d RL , set = —6, <f> = tt/2 — a. 

The case d = 0. We first argue that dub(0) = 7tt/3. The case d = is much easier since there is only 
one degree of freedom: Without loss of generality we can assume a — 0. It is easy to verify that for any /3 
there is a CCC-path of length at most 77r/3. For /3 — ir, no Dubins path has length shorter than 77r/3, 
and so dub(0) = 77r/3 (see Figurej7|. In the rest of this paper we can therefore mostly assume d > 0, and 
avoid some degeneracies. 

The three cases. We subdivide the pairs of orientations (a,/?) G T (and hence in the triangle A) into 
three cases. We define: 

A = {(a, (3) G T I d^(a,p) < 2 and d RL (a,/3) < 2}, 
B = {(a,/3) G T I d LR (a,/3) < 2 and d RL (a,/3) ^ 2}, 
C = {(a,/3) GT| d hR (a,p) > 2}. 

For clarity, let us define the parts of A, B and C that lie inside the triangle A: 

A A = A n A, B A = B n A, and C A = C n A. 
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(a) SB-path 



Figure 8: RSR-paXhs in case C. 

Note that (a,/3) G A if and only if L s r\R F ^ and R s r\L F ^ 0. Since \£ s r s \ = 2 and I^fI = 2, the 
triangle-inequality implies d L < 4 and d R < 4, so both LRL- and i?LiZ-paths exist. We will concentrate 
entirely on these two path types in case A. Note that case A does not occur for d ^ 2, as we have 
8 > d 2 R + d| L = 2d 2 + 8 cos 2 8 ^ 2d 2 . 

Similarly, we have (a, ff) € B if and only if n i?F 7^ and R$ P\ Lp = As in case A, Li?L-paths 
must exist, and d RL ^ 2 implies that i?SX-paths exist as well. We will concentrate on RSL- and Li?L-paths 
in case B. We observe the following: 

Lemma 3. A point (cr, 8) G B has 8 > tt/2, and for (a, ft) g B A we have 

sC a sC 7r/2 and 7r + a</3<27r-a. (16) 

Proof, (er, 8) <E B means d 2 R < 4 and d RL ^ 4, so > d 2 R — d RL = 8dcoso'sin<7. Since sin a ^ 0, we must 
have cos 8 < 0, and thus 8 > tt/2. So (a, 8) lies in the top half of T. This top half intersects the triangle A 



in the triangle with corners (in (a, /3)-coordinates) (0,7r), (tt/2, 3tt/2), and (0, 27r), implying (16 1. □ 



Finally (a, 0) £ C if and only if Ls fl Rf — 0- This implies that LS'-ff-paths exist. We will study LSR- 
and i?<f)i?-paths in case C. 

We can now refine our Dubins-function dub(d) by defining the following three functions: 

dub j4 (d) = sup (Q /3)eA A £(d, a, (3) - d 
dub B (d) = sup (Q (3)e b a £(d, a, /3) - d 
dubc(d) = sup (Q /3)eC A £(d, a,j3)-d 

and we have 



dub(d) = max{ dub^d), dubs(d), dubc(d)}. 



3 Case C 



Case C is the easiest, and we can handle it immediately — the rest of this paper will be dedicated to the 
study of case A and case B. The arguments presented here are already in Kim's master thesis [16j . 

In case C we have L s n Rf = 0- Since Ls and Rp are open disks, they are allowed to touch, but not 
to overlap. We show that there always exists an RSR- or L^-path of length less than d + 2tt. 

Lemma 4. For (a, (3) e C A , we have £(d, a, f3) < 2n + d. 

Proof. Since Ls n Rf = 0, the starting point S = (0, 0) does not lie in Rp, and so there is a tangent £ SR 
to Rp through S that touches Rf from above. Let a SR be the angle made by £ SR and the positive x-axis 



(see Figure 8(a) ) 
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or or 

(a) (b) 

Figure 9: LSR-pa,ths in case C. 

Let us first assume that a ^ a SR , and consider the RSR-p&th from S to F. It consists of an initial 
right-turning arc ST$, a straight line segment TsT F , and a final right-turning arc T F F, where the segment 
is tangent to Rs and Rp at the points Ts and Tp. (When a = a SR we have S — Ts ) See Figure [8] 

Let t be the vector t = TsT F , and let 5" = 5 + t. Since t = rsrp, we have Rp = Rs + t, and so 
5" lies on Rp. We claim that S' lies on the clockwise arc FTp. Indeed, any point (a, (3) £ A satisfies 
a min{/3, 2tt — /?}, which implies that t has a positive y-component. 

It follows that the length of the RSR-p&th is \ST S \ + \T S T F \ + \T F F\ = \SS'\ + \S'F\. By the 
triangle-inequality, \SS'\ sC \SF\ + \FS'\, and so the length of the i?Si?-path is at most \SF\ + 2tt = d + 2tt. 

Consider now the case where a < a SR . We show that the LSR-pa.th from S to F has length 
at most d + 2tt. The LSR-p&th. consists of an initial left-turning arc STs, a straight line segment 
TsTp and a final rig ht-tu rning arc TpF, where the segment is tangent to Ls and Rp at points Ts 



and T F . See Figure [9(a)| Here, it suffices to observe that \ST S \ + \T S T F \ < \SF\ + \FT F \, and so 
\ST S \ + \T S T F \ + \f F F\ < \SF\ + 2ir. □ 

It turns out that the bound in Lemma [4] is optimal, and we obtain: 
Lemma 5. For any d > we have dubc(rf) = 2tt. 

Proof. Lemma [4] implies that dubc(rf) ^ 27r, so it remains to provide a matching lower bound. We 
will show that i[d, ir, n) = 2n + d, and since (ir,n) 6 C A , this proves the claim. Consider a shortest 
bounded-curvature path G from (0,0, tt) to (<i, 0,7r). This path must intersect the line x = in a point p 
and the line x = d in a point q. The distance is at least d. If the path from S to p intersects Ls U i?s, 
then Ahn et al. [TBI Factl] showed that it has length at least tt. Otherwise the path avoids Ls U Rs and 
hence must have length at least tt. The same argument applies to the path from q to F, and so the total 
length of G is at least d + 2tt. □ 

Since dub(d) ^ dubc(rf), this establishes a lower bound for the Dubins cost function. 
4 Regions of the square T for < d < 2 

We will now describe the regions A, B, and C of the square T geometrically. For our purposes it will be 
sufficient to do this when < d < 2, so we assume this throughout this section. 
We define the angle 

a* = arcsin -, (17) 
and observe that for (a, j3) = (a*, 2ir — a*) as well as for (a, (3) — (tt — a* , tt + a*) we have Rs = Rp 



(Figure 10 ) 



We now argue that there is a curve (a, <5 LR (cr))o^ C r^7r in T that connects the two points (0, a*) and 
(tt, a*), lies strictly between S = a* and S — tt/2 except for its endpoints, and such that e? LR = 2 on the 



curve, g? lr < 2 between the curve and the line S = tt/2, and d LR > 2 below the curve. Recall from (14) 
that c?l R = d 2 + 4<i cos 5 sin er + 4 cos 2 5. For S — a* , we have cos 2 6 = 1— sin 2 5 = 1 — d 2 /A, and so 
d 2 TI = 4 + 4dcosa* sin a. This is equal to 4 for a £ {0, tt}, and otherwise larger than 4. For 6 = tt/2, we 
have d 2 R = d 2 < 4. Finally, for S £ (0, tt/2), we have -§gd 2 R = — 4d sin 8 sin a — 8 cos S sin S < 0, and so 
5 i y <i LR is a decreasing function for a £ (0,tt), proving the claim. 



9 



(a) a = a* , (5 = 2tt — a* 



(b) a = it — a" , = n + a* 



Figure 10: Configurations when R$ and Rf coincide. 



Consider now d Rlj = d 2 — Ad cos S sin a + 4cos 2 <5 by (15 1. For a* ^ 5 ^ 7r/2, we have cos 2 5 



1 — sin 5 ^ 1 — d /4, and so d RL ^ 4, with equality only for cr € {0, 7r} and (5 = a*. Let 

• d 

a = arcsm — , 
4 

and consider the interval ^ cr cr*. On this interval, we have ci RL 4 for 5 = (with equality only 
for (j = a*), d 2 L ^ 4 for (5 = a* (with equality only for a = 0), and d RL ^ 4 for <5 = 7r. Since for fixed cr, 
<i 2 L = 4 is a quadratic polynomial in cos S, it has at most two roots in ^ S 7r, and thus there must be 
a unique value ^ <5 RL (cr) ^ a* on the interval ^ cr ^ cr* where d RL = 2, and we have ci RL > 2 below 
the curve (5 R l(ct). The same argument applies to the interval tt — cr* ^ cr ^ 7r. 

In the remaining region cr* < a < it — cr*, ^ <5 ^ a*, we have d 2 L < 4. Indeed, in this region we 
have sin cr > cZ/4, and so d RL < -D(cr, <5), where £>(cr, 5) = d 2 — d 2 cos 5 + 4cos 2 (5. We have J^D(a, 5) = 
d 2 sine* — 8 cos S sine? = sin<5(c? 2 — 8cos<5). Since sin 2 8 ^ d 2 /4, we have cos 2 6 ^ 1 — c? 2 /4, which gives 
64 cos 2 5 ^ 64 — d 2 /16 > d 4 since d < 2. This means that <i 2 ^ 8 cos S, which implies that S -D(cr, <5) is 
decreasing. Since D(a, 0) = 4, the claim follows. 

We now exploit the fact that c? RL (cr, 5) = d LR (a, it — 8) to obtain our desired subdivision. See Figure [Tl] 

• C\ is the region <5 ^ (5 LR (cr). Inside this region we have d LR 2. 

• C2 is the region 7r — a* ^ a ^ tt , it — i5 RL (cr) ^ 5 ^ 7r. Here we have c? LR ^ 2 and d RL ^ 2. 

• C3 is the region ^ cr ^ cr*, 7r — <5 RL (cr) ^ S ^ tt. Here we have (i LR ^ 2 and d Rlj ^ 2. 

• A is the region <5 LR (cr) < S < tt — <5 LR (cr). In this region we have d LR < 2 and d RL < 2. 

• Finally, B is the remaining region, where tt — 5 LR (cr) ^ 8, but excluding C2 U C3. In this region we 
have d LR < 2 and d RL ^ 2 . 

It is clear from the description above that the five regions A, B, C%, C2, and C3 are cr-monotone, 
meaning that a line parallel to the <5-axis intersects each region in a single interval. We will also need 
that the region B A is monotone with respect to the a-direction. 

Lemma 6. For < d < 2, there are two continuous functions a <— > (3 RL (a) and a 1— > /3 LR (a) defined on 
the interval [0,a*] such f3 LR (a*) = (3 RL (a*) = 2tt — a*, and for ^ a < a* we have 

• a + tt < f3 RL (a) < 2tt — a* < f3 LR (a) < 2tt — a; 

• d RL < 2 for (3 < f} RL {a), d RL = 2 for (3 = j3 RL (a), and d RL > 2 for (3 > (3 RL (a); 

• d LR < 2 for f3 < f3 LR (a), d LR = 2 for f3 = f3 LR (a), and d LR > 2 for (3 > f3 LR {a); 

The function j3 LR is a monotonically decreasing function of a, and we have 



B 



{(a,f3) I a e [0,a*), /3 RL (a) ^ [3 < f3 LR {a)}. 



Proof. Let us fix an a € [0, a*), so the points £s an d ^5 are fixed. While (3 ranges over [0, 27r], the point 
£p makes a full circle around F. This means that the distance <i RL is strictly increasing for half a period, 
and strictly decreasing for the other half period. This implies that in the range a + TT^f3^2TT — a there 
is at most one extremum of rf RL . The same argument shows that <i LR has at most one extremum in the 
range. 
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Figure 11: The subdivision of T and A for d = 1 



a) that S 



(0, 0) lies on the boundary 
Consider now j3 = a + tx, so 8 — tt/2. By (14) 
2tt — a, so that a — n. Since a < a* 



Consider first j3 = 2tt — a*. We observe from Figure 10 
of Rf, and so e? LR < 2. Also, since a < a* , d RL > 2. 
and (15), we have d LR = d RL = d < 2. Finally, consider j3 
we have 8 > tt — a* , so sin<5 < d/2, so cos 2 8 — 1 — sin 2 8 > 1 — gP/4. By (14) and (15), we then have 
d 2 R = d 2 L = d 2 + 4 cos 2 <5 > 4, and so d LR = d RL > 2. 

Since d LR = d RL < 2 for /3 = a + 7r, d LR = d RL > 2 for j3 = 2n — a, and both functions have only one 
cxtrcmum in this range, both functions must assume the value two exactly once in this range, at values 
/3 LR (a) and /3 rl (q). For j3 = 2it — a* we have d LR < 2 and d RIj > 2, so we have /3 RL (a) < 2n — a* < /? LR (a). 
The two functions are clearly continuous, and since we have d LR = c? RL = 2 for (a, (3) — (a*, 2tt — a*) (see 
Figure [l0]), we have /3 LH (o;*) = ^ RL (a*) = 27r - a*. 



Pick any point (a,/3 LR (a)). This is a configuration where L$ and Rf are touching. If we now increase 
a, the point Is rotates left around S, and so the distance d LR increases (at least locally). This implies 
that /3 LR (a) is a decreasing function of a. 

By Lemma|3j B A is contained in the triangle 0^a^7r/2, a + 7r^/?^27r — a. Consider Figure 10 'a). 
We first observe that for a* < a ^ tt/2, the point F — (d,0) is contained in the interior of R$, and so 
d RL < 2. It follows that for (a, j3) £ B A we must have a £ [0, a*), and the expression for B A follows from 
the above. □ 



5 Explicit expressions for the length of LRL- and i?Li?-paths 

In this section we develop explicit formulas for the length of LRL- and RLR-p&ths. 

We start by a change of perspective, and consider all configurations where c! L is fixed. We choose a 



coordinate system where the line £s^f is horizontal, and is lies to the left of £f, see Figure 12(a) We 
have drawn the two unit-radius disks Af L and 7V L tangent to Ls and Lp- The points of tangency are S% 
and F£ for M L , and S\ and for iV L . Dubins [TT] showed that the length of the middle circular arc of 
a CCC-path is larger than %, and so it lies on M L . 

So any LRL-paXh first follows a leftwards arc on Ls, then switches to Af L at S%, follows the rightwards 
arc on M L until it reaches F^ , and finally follows a leftwards arc on Lp. We note that the middle arc 
on M L does not depend on the specific endpoints S and F, it is determined entirely by and F% , and 
therefore by d L . Let fi L denote half the length of the middle circular arc S^F^- We have tt/2 < /i^ ^ tt, 
and 4sin(7r — (U L ) = rf L , so that we have 

/i L = n — 4 arcsin ^ . (18) 
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(a) LiJL-paths 



(b) RLi?-paths 



Figure 12: Locations of S and F on their disks 



The same considerations apply to i?Li?-paths, see Figure [12(b) We define /x R as half the length of the 
middle circular arc of the RLR-paih, and obtain 



= 7r — 4 arcsin ^ . (19) 



We note the following angles: 



£S\l s S% = /-F^lpF\ = 2/i L - 7T, (20) 
ZSfr s S% = ZF 2 R r F Ff = 2^ R - tt. (21) 

It is now important to understand the possible locations of the endpoints S and F along the disks. 

Lemma 7. We have 

• S lies on the counter-clockwise arc S^S^ of Ls if and only if d RL < 2; 

• F lies on the clockwise arc F^Ff of Rf if and only if d RL < 2; 

• S lies on the clockwise arc S'fS'.f of Rs if and only if d LR < 2; 

• F lies on the counter-clockwise arc F 2 F^ of Lp if and only if d LR < 2. 

If d RL 2 and in addition a tt/2, then S lies on the counter-clockwise arc SqSi of Ls- 

Proof. Consider the position of Rs as S moves once around the fixed circle Ls- The center rs describes 
a circle of radius two around is- When S — , we have Rs — N L , when S — S 2 , we have Rs = M L . If 
d RL < 2, we have Rs H Lp ^ 0, and so S must lie on the counter-clockwise arc SiS 2 . If d RL ^ 2, we have 
i?s H Lp = 0, and S* must lie on the complementary arc S 2 S\. The same argument applies to i?Li?-paths 
to determine the location of F. We argue similarly for g? lr < 2 and d LR ^ 2. 

We observe next that the starting orientation at S, which is the forward tangent to Ls at S, and the 
vector SF must make an angle of a. Since F £ Lp, this is impossible for a tt/2 and S on the long 
counter-clockwise arc S^Sq, and so d RL ^ 2 with a 7r/2 implies that S 1 € SoS^. □ 

Lemma allows us to add another case to Lemma [TJ 

Lemma 8. Let d\ < d 2 , and (a, (3) £ B A (where B A is defined for d\). If there is an LRL-path of 
length I from (0, 0, a) to (d\, 0, f3), then there is a path of length £ + (d 2 — d\) from (0, 0, a) to (d 2 , 0, /3). 

Proof. We only have to prove that the LiiL-path has a horizontal tangent oriented in the positive 
x-direction. Assume this is not the case, so there is no point on the path were the orientation is or 2tt. 
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The path starts at orientation a, the orientation increases to a + 71, decreases to a + 71 — 2/i L , and 
increases again to f3 — a + 71 — 2/i L + 72, without ever leaving the open range (0, 2tt). But by Lemma [3j 
(a,/3) € F A implies d RL ^ 2 and ^ a ^ 7r/2 and thus, by Lemma[7]and (20), we have 2/i L — tt ^ 71 ^ /i L 



and ^ 72 ^5 2/i L — tt. This implies that 71+72 ^ 2/x L , and so (3 ^ a, a contradiction to (a, (3) £ B A . □ 
Lemma 9. For any a, 5, we have 

LRL(ct, 5) = Afi L + 25 (mod 2tt), 
RLR(ct, 5) = 4yu„ — 26 (mod 27r). 

For (a, 6) £ A, we have 

lrl(ct, <5) = 4/i L + 25 - 27r, 
RLr(ct, <5) = 4yu„ - 25. 

For (a, 6) £ B, we have 

lrl(ct, 5) = 4/i L + 25 ~ 27r, 
rlr(ct, <5) = 4/i„ - 25 + 27T. 

Proof. An LFFpath consists of an initial left-turning arc of length 71 on Lg, a right-turning arc of 
length 2/i L on the middle disk, and a final left-turning arc of length 72 on Lp . This means that the total 
change in orientation is 71 — 2/i L + 72 . On the other hand, since the initial orientation is a and the final 
orientation is /?, this must be equal, up to multiples of 2tt, to (3 — a = 26. It follows that 

lrl = 71 + 72 + 2 ) u L = 71 - 2/i L + 72 + 4/x L = 25 + 4jU L . 

For FLF-paths, we can similarly observe that —71 + 2/i R — 72 = 25 (mod 2tt) (here, 71 and 72 are the 
right-turning arcs) and obtain 

RLR = 71 + 72 + 2fi„ = 4/i R - (-71 + 2/i R - 72) = 4/i R - 25. 

Let us now assume that (a, 5) £ A. We have 2a* < 26 < 2it — 2a*. On the other hand, 71 + 
72 - 2/i L ^ -2^ L ^ -27T. By Lemma [7J S £ S^Si; and F £ F 2 L F 1 L , and so we can extend the LRL- 
path to a complete clockwise loop as in Figure p(b)| The loop uses additional left-turns (1 and £2, 
and an additional right-turn of length 2/i L . The total turning angle of a clockwise loop is — 2tt, and 
thus 71 + 72 + 5i + 62 — 4/i L = — 2tt. Since 2/i L ^ 2tt this implies that 71 + 72 — 2/i L ^ 0. From 
-2tt < 71 + 72 - 2^ L < and < 2a* < 25 < 2tt - 2a* ^ 2tt, we conclude that 71 + 72 - 2/i L = 25 
(mod 2ir) implies 71 + 72 — 2fi L = 26 — 2ir. This shows that lrl = 4/i L +26 — 2ir. 

For FFR-paths, we could argue analogously, or we can simply observe that 

rlr(<t, 6) = lrl(7t - a-, tt - 6) = 4^ L (7r - cr, tt - 6) + 2(tt - 5) - 2ir = 4/i R (cr, 5) - 25. 
Assume now that (a, 5) £ B. By Lemma |3j we have tt < 25 ^ 2tt. By Lemma [7] and (20) we have 

2/i L — TT ^ 7! < 2-7T 

^ 72 < 2/i L - 7T 

Together these imply —tt ^5 71 + 72 — 2/i L < tt. Since 71 + 72 — 2/i L = 25 (mod 2tt), we must have 
7 + 72 - 2/i L = 25 - 2-7T. It follows that lrl(ct, 5) = 4,u L + 25 - 2tt. 

We turn to the FLF-path in case B. By Lemma [7] and (21) we have (here, 71 and 72 are the 
right-turning arcs) 

71 2/i R - TT 
2^i R - tt 72 < 2tt 

Together we have — tt ^ 71 + 72 — 2/i R ^ tt. Since —71 — 72 + 2/i R = 25 (mod 27r), we must have 
—71 + 2fi R — 72 = 25 - 27r, which shows that RLR = 71 + 2fi R + j 2 — 4^r — (2(5 - 2tt) = Afx B —25 + 2tt. □ 
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(a) ^ 71 , 72 ^ — 1" (b) Completing to a loop 



Figure 13: Proof of Lemma [9] for case A 
6 CCC-paths for d < 2 and case A 



When d < 2, both LRL- and RLR-pa.ths exist for any (a, 8) € V. In this section we analyze the length of 
these two CCC-paths for < d < 2. 

We define three functions L, R, and c on T: 



L(cr, 5) = Afi L (a, S)+28-2tt 

r(ct, 5) = 4/i R (cr, (5) - 25 

c(cr, 5) = min{L(cr, (5), r(ct, 5)}. 



(22) 
(23) 
(24) 



While these functions are defined and continuous everywhere on T, we have shown in Lemma [9] only that 
LRL(ct, 8) — l(ct, 8) for (a, 8) e AlJ B, and RLR(cr, S) = r(ct, 5) for (a, 8) G A. It follows that c(cr, <5) is the 
length of the shortest CCC-path for (a, 5) e A. 



We obtain the derivatives of L and R using /i L 



- arcsin ^ , /i R = 7r — arcsin % , and Equations ( 12 )- 



dh 
OS 
dh 
da 
<9r 
88 
Or 
~0a 



{a, 5) 
(a, 5) 
(a, 8) 
(a, 8) 



-4 cos S sin 8 + 2d cos 8 cos c 



Dr. 



-2G?sin<5sin<7 



D, 



-4 cos (5 sin 8 — 2d cos 8 cos er 



D R 



2cZ sin 8 sin (j 

d7~ 



(25) 
(26) 
(27) 
(28) 



The derivatives are not defined when D L = or D n = 0. D L = occurs when <i L = or <i L = 4, and 
D R = occurs when d R = or d R = 4. 

Recall that a* — arcsin |. We define the following rectangle S C F: 



$J a < 7r and 



< 5 € 7T - 



Note that icH (see Figure 11(a) ). In the interior of S, we have sin 8 > d/2, which implies cos£ 2 < 
1 — d 2 /4. So for a* < 8 < tt/2, we have d RL < 2 by (15), and for 7r/2 ^ 5 < ir — a*, we have 
d LR < 2 by (14). By the triangle inequality it follows that d L < 4 and d R < 4. Also, by (12) d\ = 
(d — 2sin(5coscr) 2 + 4sin<5(— coser + 1). Since sin<5 > 0, d L = can occur only when a = and 
8 € {a*, 7r — a*}, which occurs at a corner of 3. Similar arguments hold for the case where d R = 0. Thus, 
in the interior of S, D L ^ and D R 7^ 0. 



14 



Lemma 10. For (a, 8) € H, the function 

• o~ i y l(ct, 5) is decreasing, 

• a i — y r(ct, <5) is increasing, 

• (5 i ^ l(<t, <5) is increasing, 



8 i— > r(ct, (5) is decreasing. 



Proof. Consider Equations (26) and (28). Since sin 5 > in S and sin a ^ with equality only for 
o~ € {0, 7r}, we have < and ^ > in the interior of S. 

It remains to discuss the two functions of 5. We show that §§(cr, 5) > 0. If Z = 2cos5(— 2 sin 5 + 
d cos cr) ^ 0, this is true. Let us thus assume that Z < 0. Since sin 5 ^ d/2 implies that — 2sin5 + dcoscr ^ 
—2(d/2) + dcosa = d(— 1 + coser) ^ 0, we have cos5 > and so — 2sin5 + dcoscr < 0. sin5 ^ d/2 
also implies that d 2 ^ 4 sin 2 5, so by ( 12 ) d 2 — d 2 — Ad sin 5 cos cr + 4 sin 2 5 ^ 8 sin 2 5 — 4d sin (5 cos a = 
A sin 8(2 sin 5 — dcoscr). Since d 2 ^ 0, we have 

d\ ^ 16 sin 2 5(2 sin 8 — dcosa) 2 . 

On the other hand, we have 

d 2 = d 2 - Ad sin 5 cos a + A sin 2 5 

^ d 2 cos 2 a — Ad sin 5 cos cr + A sin 2 (5 = (2 sin 5 — d cos a) 2 . 



(29) 



(30) 



Now we want to show that 2Z? L ^ — Z, which will imply §j(cr, 5) 0. Let us compare the squared terms: 



AD 2 - Z 2 



AD 2 - A cos 2 5(-2 sin 5 + d cos a) 2 

dV 



= 4 U 1- 



16 



— cos 6(2 sin 5 — d cos a) 



^ 4(d 2 - sin 2 5(2 sin 5 - dcoscr) 2 - cos 2 5(2 sin 8 - dcosa) 2 ) 
= 4(dj? - (2 sin 5- d cos cr) 2 ) ^ (by (j30j|). 



(by (29)) 



One of the inequalities in above formula is a strict inequality: if ( 30 ) is an equality, then d 2 = (2 sin 5 — 
d cos cr) 2 , which means that cos 2 a = 1. This implies that d 2 < 4sin T 5 since we have — 2 sin 5 + d cos cr < 0, 
so (29) is strict. In the case where |29| ) is an equality, we can argue similarly that (30) is strict. 

Similarly, we prove that ||(cr,5) < 0, since sin5 ^ d/2 again implies that d 2 = d 2 + 4dsin5cos a + 
4 sin 2 6 ^ 4 sin 5(2 sin 8 + dcosa), so 



and we have 



d R < 16 sin 2 5(2 sin 5 + d cos cr) 2 , 



d 2 ^ d 2 cos 2 cr + 4d sin 5 cos cr + 4 sin 2 5 = (2 sin 5 + d cos cr) 2 



□ 



We also need to show that the length of Li?L-paths is monotone in a and f3, at least for the cases of 
interest to us. 

Lemma 11. For (a,j3) 6 B A , the function a t— > LRL(a, /3) is decreasing. 

Proof. By Lemma|9]we have lrl(q;,/3) = 4/i L + 25 — 2n in £? A . Since 25 = ft — a, /« L 
d 2 = (d - sin /3 + sin a) 2 + (cos /3 - cos a) 2 from |6} , and 8 a ™ n x = ,} n , we have 



7r — arcsm ■ 



d lrl (d — sin (3 + sin a) cos a + (cos (3 — cos a) sin a 

-doT {a ^ ) = Dl 



We also have 



(d — sin j3 + sin a) cos a + (cos ft — cos a) sin a =d cos a — sin /3 cos a + cos /3 sin a 

=d cos a — sin(/3 — a) ^ 0. 

The last inequality holds since a it/ 2 and j3 — a ^ 7r in _B A by Lemma [5} It follows that 
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Lemma 12. For (a, j3) £ A A U B A , the function f3 i ^ LRL(a,/3) is increasing. Moreover, if (a, /3) £ B A 
and then ^p(a,/3) > 1. 

Proof. We have 

<9lrl (d — sin/3 + sin a) cos/3 + (cos/3 — cos a) sin/3 

= Dl + 

Note that the derivative is not defined if -Dl =0, but as we noted before, that can happen only on the 
boundary of the region. 

We first prove the second statement. /3 ^ ^ implies that sin/3 ^ and cos/3 ^ 0. By ( [l6| we have 
7r + a^/3^27r — a, which implies that cos/3 ^ cos a. It follows that the first term of (3) is 

nonnegative, proving that ^p(a, /3) ^ 1. 

Let U = d — sin /3 + sin a, and V = cos j3 — cos a (note that d 2 = C/ 2 + V 2 ) . If U cos (3 + V sin /3 ^ 0, 
we have i =^|p(a,/3) > 0, and /3 H> LRL(a, /3) is increasing. Now let us assume that [/cos/3 + Vsinft 

is negative, and let us compare the squared terms of d gp L (we only consider the numerator since the 
denominator is always positive). 

/ d 2 \ (U 2 + V 2 ) 2 

(U cos(3 + V sin /3) 2 - d 2 1 - -i =(f/ 2 cos 2 /3 + F 2 sin 2 /3 + 217V cos /3 sin /3) - (U 2 + V 2 ) + i - 

\ 16/ 16 

(•7-7-2 , T/2\2 

=(-U 2 sin 2 /3 - I/ 2 cos 2 /3 + 2UV cos /3 sin /3) + - '- 

16 

= — (—16(17 sin ^ - y cos /3) 2 + (C/ 2 + y 2 ) 2 ) . 
16 

For (a, /3) e A A U B A , we have d LR < 2. From ([8]), d 2 R = (d + sin /3 + sin a) 2 + (cos /3 + cos a) 2 < 4. 
By substituting d + sin a = U + sin /3 and cos a = cos /3 — V , we have 

^(f7 + 2sin/3) 2 + (2cos / 3-V) 2 <4 

<S> C/ 2 + V 2 + 4(Z7 sin /3 - V cos 0) + 4(sin 2 /3 + cos 2 /3) < 4 

^ (C/ 2 + V r2 ) + 4(C/sin/3-ycos/3) < 0. 

Since U 2 + V 2 = dl ^ 0, this implies that U sin/3 — ycos/3 < 0, and further the squared term 
16{U sin /3-V cos f3) 2 is greater than (U 2 + V 2 ) 2 , which implies that (U 2 + V 2 ) 2 - 16(U sin /?- V cos /3) 2 < 0, 
completing the proof. □ 



Lemma 13. TTie function c(o~,5) has no local extremum in the interior of 3 except at {tt/2, 7r/2). 



Proof. By Lemma [lOj neither L nor R has a local extremum in the interior of 2, so any local extremum 
of c(a,S) must be a point in the set A of points (a, 5) with l(<t,S) = r(ct, £). By Lemma 10 A is a 
(5-monotone curve. Since l(ct, <5) = r(tt — a, tt — 8), the curve A passes through the point (tt/2, tt/2). By 



Lemma 10 this implies that l(o~,5) < r(ct, 8) for the quadrant tt/2 ^ a ^ tt, oj* ^ 8 ^ 7r/2, and that 
r(ct, <5) < l(<t, 5) for the quadrant ^ a ^ 7r/2, 7r/2 ^ r5 ^ 7r — a* except at the point (tt/2, tt/2). By 
point symmetry, we can restrict our attention to the range tt/2 < a < tt, tt/2 < 8 < tt — a* . 

Assume for a contradiction that (<r, 8) £ A is a local extremum of L, restricted to A. This implies that 
the gradient Vl(ct, 8) and the normal of A in (o~,6) are linearly dependent, by the method of Lagrange 
Multipliers. The normal of A is the gradient of h(a, 8) — r(ct, 8), so Vl(ct, 8) and Vr(<t, 8) must be linearly 
dependent. Note that Vl(ct, 8) = (f^(c, 8), §§(c, 8)), and Vr = (§f (o - , 8), §§(ct, <£)), which we can obtain 



from Equations (25 )— ( 28 ) 



For the two vectors to be linearly dependent, we would have to have 

<9r , 



D^io-,8) 



D R ~-(a,8) 



0, 



which means 



2L> L - 2D R - 8 cos 8 sin 8 = 0. 



16 



7T 




Figure 14: The curve A where h(a, S) = r(<t, 5) in rectangle S (shaded region), shown for d = 1. 

In the range under consideration, — 8 cos S sin S > 0. We will show that D L > D R , a contradiction. We 
have 



16(Dl 



Dl) = dt(W-dt)-dl(W-di) 

= di -4 + i6(dl - 4) 

= (dt-dl)(16-(dt + dl)) 

= — 8d cos (7 sin (5(16 — (2d 2 + 8 sin 2 5)). 



Since coscr < and d < 2, the expression is positive. 



□ 



Since c(a,S) is continuous, it assumes its maximum on H. By Lemma |13| this must happen either 
at (a, S) — (7r/2, 7r/2), or on the boundary of S, at a point (er A , S A ) where l(<t a , 5 a ) — r(ct a , <5 a ). As we 



observed in the proof (see Figure 14), this must happen either on the vertical side a — tt,tt /2 < 5 ^ tt ~ a* , 
or on the horizontal side 6 — ir — a* , 7r/2 < a ^ 7r. The point (<t a , 8 A ) is unique if we require a A + S A > 7r, 
and there is a symmetric point (7r — er A , n — S A ). Let us define the function a(cQ for ^ d < 2 as 

A(d) = lrl(ct a , S a ) = rlr(ct a , S a ). (31) 
There is an important breakpoint at d — \/2: 



Lemma 14. The maximum A(d) occurs with o~ A 
V2^d<2. 



7r when ^ d ^ and with S A = ir — a* when 



Proof. We evaluate l(7t, it — a*) and r(7t, n — a*). Using ( 12 ) and ( 13 ), we have 



which implies by (22) and (23) 



<i L (7r, tt — a* 
d R (TT, tt — a* 



l(7t, tt — a*) 
r(7t, tt — a*) 



2d 
0, 



Att - 
2tt- 



6a*, 
2a*. 



Since a* = arcsin(d/2), we have r(tt,tt — a*) < l(tt 7 tt — a*) for d < \^2, equality for d = \[2, and 
R(7r,7r — a*) > l(tt,tt — a*) for d > \[2. In the first case, Lemma 10 implies that the maximum must 
occur on the vertical side at a = tt, in the last case it must occur on the horizontal side at 5 — tt — a*. 
For d = \/2 the maximum occurs at the corner (aA, Sa) = (tt,tt — a*). □ 

Lemma 15. On the interval ^ d ^ \[2, the function 
• d i y A(d) is monotonically increasing, 
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• d i y A(d) — d is monotonically decreasing. 

Proof. We will show below that the two functions d i— > A(d) and d n- A(d) — d have no extremum on the 
interval (0, v2). This will imply the claim if we observe that 

A(0) = 7tt/3, 

A(v / 2) = 5tt/2 > A(0) 

A(V2) - V2 = 5tt/2 - V2 < A(0) - 0. 

Again we will employ Lagrange multipliers. Let us first give the necessary derivatives. Setting a = tt 
and 7r/2 < 5 ^ tt — a* (by Lemma 14), we have: 

d L = d + 2 sin <5 
d R = 2 sin 5 — d 



by (12| 



by ([13J) 

L(d, 7r, S) = 2tt - 4 arcsin(d L /4) + 25 
R(d, tt, 5) = 47r - 4 arcsin(d R /4) - 2<5 



by ([221 
by pi 



Let us introduce F L = v/l — (d L /4) 2 and F R = y 7 ! — (d R /4) 2 to obtain: 



_5 
9d 

dd 



L(d, 7r. (5) 
R(d, 7r. 6) 



1 

'F 



L(d, TT, S) 



1 



9 



R(d, 7T, (5) 



2- 
-2 



2 cos <5 

F L 
2 cos 5 



dd " /•',. 

We consider first the function d H > A(d). An extremum of A(d) is an extremum of the two-parameter 
function (d, <5) >-> h(d, tt, S) under the restriction h(d, tt, 6) = R(d, tt, S). Such an extremum would have to 
satisfy the condition VL(d, tt, 5) = AVR(d, tt, S). For this to hold, 

2 - 2cos<5/F L 



F L -2-2cos£/F R ' 



which implies 



F n - F, 



-2 cos S. 



(32) 



But in the following we show that this is impossible. Since we have 



. d L . d R 

:csm arcsm — 

4 4 

for the condition h(d, tt, 5) = R(d, 7r, <5) to hold, we need 



>4 



dR 
4 



2d, 



27r - 45 < -2d S ^ - + - sin (5 < cos - 
2 2 2 

^ • d 

— cos o ^ sm -. 

2 



(33) 
(34) 



We also claim that F R + F L > 1 ; indeed since d R < 2, F R > V3/2, and since d L < 2 + v/2, F L > 1 - V3/2. 
Since F 2 - F 2 - (F R - F L ) = (F R - F L )(F R + F L - 1) ^ and F R = F L only if d = or 5 = tt/2, this 
implies that 



F R -F L <F 2 -F 2 = ^ sin ,5. 



(35) 



Now, 



F r! 



F L < - sinJ 
d d 

^ 77 COS - 

2 2 

„ . d 

< 2 sm - 



by (351 



by (331 



^ —2 cos 6 



by (34|. 
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The third inequality holds since 

d fd d d\ 1 d d d d 

— - cos 2 sm - = - cos sin cos - CO, 

dd\2 2 2) 2 2 4 2 2 ' 

which implies that D(d) = (d/2) cos(g?/2) — 2sin(d/2) is a decreasing function, and finally note that 
D(0) = 0. 

Consider next the function d M- A(d) — d. An extremum of A(d) — d is an extremum of the two- 
parameter function (d, S) H ► L(<i, n,S) — d under the restriction l(cZ, 7r, o") = R(d, 7r, (5). Such an extremum 
would have to satisfy the condition XV(h(d, n, 5) — d) = V(l(c£, tt, 6) — R.(d, tt, 5)) , or 

1 ^ ^ 2 cos 5 \ / 1 1 ^ 2 cos <5 . 2 cos <5 



F, ' F L y V F L F R ' F L F R 

The two components give us the following conditions on A: 

1/F R - 1 



A = 1 
A = 1 



1/F L + 1 

2+ (2cos<5)/F R 

2 - (2cosc5)/F L 



This is equivalent to 

/ 2 cos 6 \ / 1 \ / 2 cos S \ / 1 
Multiplying out and rearranging the terms gives 

Since cos 5 < and F R > F L , the left-hand side is negative. We will now show that the right-hand side is 
positive, a contradiction, and so d i— > A(d) — d cannot have a local extremum. 
It is enough to show that cos 2 8 ^ (F R F L ) 2 : 

(^•-.^.-(i-i+S + fiW!)-.^. 

16 16 2 

-2d 2 + 8 sin 2 5 (d L d R ) 2 



> 0. 



16 16 2 

The last inequality holds since sini5 ^ d/2. 

It remains to decide whether c(d, tt/2, 7r/2) or A(d) is larger. 
Lemma 16. For < d ^ a/2, we ftawe max^s)es c(d, a, 5) — A(d). 
Proof. For < d s? \/2 and (cr, <5) = (tt/2, tt/2), we have 

dl=d 2 + 4 by @ 
/i L = 7T — arcsin(\/ d 2 +4/4) 
L(d, tt/2, tt/2) = 4/i L + 25 - 2tt = 3tt - 4 arcsin(Vd 2 + 4/4). 
Since arcsin is an increasing function, d i— > L(d, 7r/2, 7t/2) is a decreasing function. We therefore have 

L(d, tt/2, tt/2) < l(0, tt/2, tt/2) = 7tt/3. 



□ 



On the other hand, by Lemma 15 the function d H> A(d) is increasing, and so A(d) ^ a(0) = 7tt/3 > 
L(d, tt/2, tt/2). □ 
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We now have all the tools to discuss case A. In case A, which occurs only for d < 2, we have d LR < 2 
and d RL < 2. We will now justify that it suffices to study CCC-paths in this case, as no other path type 
can be shorter. Since LSR- and i?5X-paths do not exist, it is enough to show the following lemma: 



Lemma 17. For (cr, 8) € A, we have 



lrl(ct, 8) < lsl(ct, 8) 
rlr(ct, 6) < RSr(<t, 8). 



Proof. Let 71 and 72 be the length of the left-turning arcs of an LRL-uath. By Lemma [7J the endpo 



nts S 
The 

LSL-pa,th turns left on L$ until Sq, goes along the tangent to Fq, then turns left on Lp until it reaches F. 

g? l , we have 



and F lie on the counterclockwise arcs S^S^ of Lg and F^Ff of Lp. Consider now Figure 12(a^ 
LSL-pa,th turns left on L$ until t 
Since ZS%£ S S% = ^ and |S£F L 

lsl - lrl = c? L + 2(27r - ,u L ) - 2^i L = d L + 4(7r - ,u L ) > 

since ^ L ^ tt. The analogous argument shows that RLR ^ RSR. 



□ 



Lemma 18. For < d < \[2, dub y i(d) = A(d) — d. In other words, the maximum is realized by the 
unique point (o~ A ,8 A ) on the segment a — u, n/2 ^ 8 ^ ir — a* where h(d,a A ,S A ) — R(d,o~ A ,S A ). 



Proof. By Lemma 
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we have max( CT .5) 6 h c(d, a, 8) — a(gQ, and the maximum is assumed at the point 
mma 17 we have £(d, a, 8) = c(d,a,8) for 
this means that dubA(d) = sup^^^^A £(d, 0,8) — d = c(d, a A , 8 A ) — d 

Lemma 19. For \[2 ^ d < 2, we have dub^d) ^ max{27r, dubs(d)}. 



{o- A ,5 A ) € A A . By Lemma [l7| we have t{d,a,8) = c(d,a,S) for (a, 5) <E A A . Since (a A ,S A ) € A A C 5, 

A(d) -d. □ 



Proof. Let A be the closure of A A . Since A is compact and C is continuous, there is a point (cr, 8) G A 

where l(ct, 5) = r(ct, <5), 
it cannot occur on the 



where c(o~,8) assumes its maximum. By Lemma 13 this is necessarily a point where l(ct, 8) = r(ct, 8), 
and either (cr, 8) = (tt/2, 7r/2), or (a, 8) lies on the boundary of A. By Lemma 
vertical side of S. 
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Assume first that 8 < n/2. By Lemma 10 we must then have a < tt/2. Using Lemmas 10 and 12 
have 



we 



c{d,a,S) ^ h(d,a,8) < h(d,8,8) ^ h(d, tt/2, tt/2). 



We observed in the proof of Lemma 16 that h(d, 7r/2, 7t/2) is a decreasing function of d. For d = \[2, we 
already have l(a/2, vr/2, tt/2) — 3ir - 4 arcsin(v / 6/4) < 27r + \/2, and so c(d, a, 8) ^ d+ 2tt. The same 
argument covers the case where (cr, 8) = (tt/2,tt/2). 

It remains to consider the possibility that 8 > tt/2. In this case (cr, (5) lies on the common boundary 
of A A and B A . On this boundary we have d RL = 2, so Rs and Lf touch. Note that in this case (cr, 8) 
lies in i3 A , not in A A . It remains to observe that i?LiZ-path is identical to the -RSX-path, so we have 
RSL(cr, 8) = lrl(ct, 8). We will prove in Lemma |23| that in B A these two path types are always shortest, 
and so c(d, cr, 8) ^ dub^ (d) + d. □ 



7 i?5X-paths for d < 2 and case B 

It was proven by Goaoc et al. [13] that for any CSC-path type (that is, one of the types LSR, RSL, 
LSL, or RSR), the length of a path of this type from (0, 0, a) to (d, 0, /3) is differentiable at any point 
(a, fi) £ □ where such a path exists and both its circular arcs have non-zero length. For the case of 
RSL-pa,ths, they prove specifically that 

1 — cos 7 R (36) 
l-cos7 L , (37) 
where 7 R and 7 L are the lengths of the right-turning and the left-turning circular arc on the path. 



d_ 
da 
d_ 
88 



RSL(a, B) = 
RSL(a, P) = 
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Figure 15: When 7 L vanishes. 



We recall that case B is the situation where c£ LR (a, (3) < 2 and e? RL (a, j3) ^ 2. For < d < 2, Lemma[6] 
gives an explicit description of the region B A , using the two functions /3 RL (a) and /3 LR (a). Let us define 
two extended regions: 



B° = 
B = 



{(a,p) |0sCasCa\ /3 RL (a) < /3 sC 2tt - a}, 
{(a,/3) lOsCasCa*, /3 RL (a) < /3 < /3 LR (a)}. 



(38) 
(39) 



So .B is the closure of B A , while B° is the union B A U C£ (see Figure 11(b)) 



We now investigate where the three segments of an i?SX-path can vanish in B°: First, the ^-segment 
vanishes exactly if d RL = 2. This happens exactly on the lower boundary of B°. By Lemma [7j S lies 
on the arc SqS^ of L$ (see Figure |l2j), and so < 7 R ^ 2a* . If 7 R = 0, then we must have S — Sq and 
therefore F = F^. This is the case (a, j3) = (0, 2n). Finally, if 7 L = 0, then F must lie on the arc F£ Fq , 
and we have 2tt — /3 ^ a (see Figure [T5] ) . Equality holds only for S — , F = F^ , which is the case 
(a,/3) = (a*, 2-7T — a*). In all other cases, 2ir — (3 < a is a contradiction to (a, /3) € A. 

The above implies that RSL(a,/3) is differentiable in any point in the interior of B°. The function 
is continuous everywhere except at the two points (a, j3) = (0, 2ir) and (a, f3) — (a* , 2ir — a*). At the 
first point, the i?«SX-path degenerates to the line segment SF of length d, while the limit of RSL(a,/3) 
for (a, (3) — > (0, 2tt) is d + 2n. For the second point (a*, 2n — a*), consider Figure [To](a). At this point, 
both the straight segment and the left-turning arc vanish at the same time, and the length of the path is 
only \SF\ = 2a*. However, for (a,/3) -> (a* , 2tt - a* ) , the limit of RSL(a, /3) is \SF\+2ir = 2a* +2ir. We 
observe that this is exactly the value of RLR(a*, 2ir — a*), as the final right-turning arc of the i?Li?-path 
vanishes. 

We therefore define the following function on B°: 



RSL*(a,/3) 




for (a, (3) 
for (a, (3) 
else 



(0,2tt) 
(a*,27r 



(40) 



We have 

Lemma 20. For (a, f3) 6 B° , the function 

• a i — y RSL*(a,/3) is increasing, 

• {3 i — y RSL*(a,/3) is increasing, 

• a i— > RSL*(a, 27r — a) is increasing, 

• RSL* (a, /3) < RSL* (a* , 27T - a*) = 2a* 



2tt. 



Proof. The two derivatives ( 36 1 and ( 37 1 are defined and both positive in the interior of B° , implying 
the first two claims. 

For the third claim, we need to show that 



sC — RSL(a,/3) 
oa 



d(3 



RSL(a,/3) 



cos 7 R + cos 7l 



for (3 = 2tt — a. By Lemma [TJ S lies on the arc SqS^ of L$. If (3 = 2tt — a, then F lies on the 
counter-clockwise arc F^Fq of Lp. The two circular arcs of the i?5X-path have two components, namely, 
7a = a + C and 7 L = /3 + C = 27r — a + C, see Figure [16] This implies 
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Figure 16: An i2SX-path for j3 = 2ir — a 

hs 




Figure 17: If 7 R + 7 L ^ 2n then rsl ^ d + 2tt 
— cos 7 R + cos 7 L — — cos(a + C) + cos(a — Q = 2 sin a sin ( ^ 0. 
The second and third claim immediately imply the last one. 



□ 



If the two circular arcs of an i?5X-path add up to less than a full circle, then the path cannot be too 
long: 

Lemma 21. For (a, j3) e B° and d < 2, if j R + >y L < 2it, then RSh(d, a, /J) d + 2tt. 

Proof. Let T$ and T F denote the points of tangency of the ^-segment to Rs and L F . We observed above 
that 7 R ^ 7T. If we also have 7 L ^ 7r then d ^ |TsTp|, and the claim follows immediately. 

If 7r ^ tt/2, then we have 7 L ^ 2tt — -f R ^ 3ir/2. But then d = |£.F| ^ 2, a contradiction. It follows 



that we must have 7 R < tt/2. See Figure 17 Let S' be the point on Lp such that the counter-clockwise 



arc S'Tp on Lp has length 7 R . Let hs and /ig/ be lines through S and 5' orthogonal to the segment TsTp. 
The distance between hs and hs' is |TsTp|, and so we have \SS'\ ^ |TgTp|. By the triangle inequality, 
we have d+\FS'\ = \SF\ + \FS'\ ^ \SS'\ ^ \T S T F \, so \T S T F \ - \FS'\ ^ d. Since 2ir - \FS'\ = 7r + 7l , 
we have 

RSL(d, a, (i) = 7r + \T S T F \ + 7l = \T S T F \ +2tt~ \FS'\ ^d + 2n. □ 
Lemma 22. If a i— > RSL*(d, a, /3 LR (a)) has an extremumforO < a < a*, i/ien RSL*(rf, a, /3 LR (a)) ^ d+27r. 
Proof. Let ao be such an extremum, and consider the i?SX-path for (ao, /3 LR (ao)). If 7 R + 7 L ^ 27r, then 



by Lemma 21 we have RSL(d, a , /3 LR (a )) < d + 2tt. We can therefore assume 7 R + 7 L > 2n. The point 
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(acb /3lr(c*o)) is an extremum of the function RSL*(a,/3), under the constraint that d^ B = 4, and so there 
must be a constant A such that V RSL* (a 0) /?lr(Qo)) = XVd% R (&o, /3lr(«o))- 
Using tt/2 < a,5 ^ tt and (14) we have 



do 



4<icos5coser > 

—4 sin 6(d sin tr + 2 cos 5) . 



Using ( 36 ) and (37 1 we get 



2^" RSL* 

do 

2 ix RSL * 
do 



= ^ RSL 



d_ 
da 

d_ 
da 



RSL 



2 — cos 7 L — cos 7r ^ 



RSL* = — COS 7 L + COS 7r. 



(41) 
(42) 

(43) 
(44) 



Inequalities (41) and (43) imply that A ^ 0, so let us consider (42). By Lemma [6j we have o + S = 
/3lr(c*o) ^ /?lr(«*) — 27r — a* > 37r/2 for d < 2, and so cos 6 < cos(37r/2 — a) = — siner. It follows that 



2 cos S < — 2 sin rr < — d sin o, and so (42) is positive 



We have ^ 7 R ^ tt and we assumed that 7 R + 7 L > 2tt. It follows that 7 L > 2tt — 7 R ^ tt. Since 
7 R > 27r — 7 L we have cos7 R < cos7 L . This implies that (44) is negative. But this means that A < 0, a 
contradiction. 



□ 



Let us define the following function B(a,/?) on B°: 
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and 
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(45) 

B{d) - d. 

we have 



B(d, a, f3) = min{RSL* (d, a, /?), LRL(d, a, (3)}. 

Our goal will be to determine b(c£) = sup( a m 6B A b(g?, a,/3), and then to show that dubs(d) 
Since B(d,a,j3) is continuous, we have s(d) — max(„ )l g)gs B(d, a, /3). By Lemmas 
B(d, a, /3) ^ B(d, a, /9 LR (a)) for any (a, 0) e -B, and so 

B(d) = max _ B(d, a, /3) = max B(d, a, /3 LR (a)). 
(a,p)eB O^a^a* 

We now show that dubs(cQ = b(g?) — d by showing that the LRL-paXh or the TiSX-path is shorter 
than any other path type. Since there is no L5"iZ-path in case B, it suffices to exclude path types LSL, 
RLR, and RSR. 

Lemma 23. For (a,/3) € B A we have 

LRL(a,/3) LSL(a,/3) 

LRL(a,/3) KLS.(a,0) 

RSL(a,/3) < RSR(a, (3) 

Proof. We first compare the Li?L-path and the LSX-path. Let 71 and 72 be the two arcs on the LRL- 
path. By Lemma [7J we have 2/i L — tt 71 ^ /i L and ^ 72 ^ 2/i L — 7T. The LSX-path has length 
LSL = 71 + 72 + 2(27r - /i L ) + d L . We thus have LSL - lrl = d L + 4(ir - /i L ) > 0. 

Consider now the Li?L-path and the i?Li?-path. By |l2| and |l3]) , we have d B —d^ = 8d sin S cos a ^ 
for (<r, <5) 6 -B A , and so d R ^ c? L , implying ^ L ^ /^ R . By Lcmma|9j we have RLR — lrl = 4(/z R — /z L — <J+7r) ^ 
0. 

Finally, we compare i?5X-path and RSR-p&th. The i?SX-path consist of an initial right-turning arc 
STs, a segment T$Tp, and a final left-turning arc TpF . The RSR-p&th c ons ist of an initial right-turning 

2, by Lemma [7J the 
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arc SRi, a segment R1R2, and a final right-turning arc R2F, see Figure 

We first claim that the arc STs is common to both paths. Indee d, sinc e e? LR < 
initial arc SRi of the i?S'i?-path must have length at least tt (Figure 12(b) ), while the arc STs must be 



shorter than tt (Figure 12(a)). 
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Figure 18: RSL(a,/3) < RSR(a, /3) for (a,0) E B A with /3 close to /3 LR (a). 



Since \T F F\ = 2tt - |FT F | and \T s Ri\ = 2n - \RiT s \, where FT F is a left-turning arc on the disk L F 
and R\Tg is a right-turning arc on disk R$, we have 

RSL - RSR = (|7yZ>| - |FT F |) - 



LRiT s U T S T F \ - \R X R 2 U R 2 F U FT F | 



The path R1R2 U i? 2 / U FT F is a path connecting i?i with T F while avoiding the interior of Rs- However, 
the path R\Ts U TsT F is clearly the shortest path of this kind, and so RSL — RSR ^ 0. □ 

It remains to understand the function B(d). By definition, there is an a B E [0, a*] such that 
B(d) = B(d,a B ,/3 LR (a B )). 

Lemma 24. The function a 1— > LRL(a, f3 LR (oij) — RSL(a, /3 ii? (a)) is monotonically decreasing on [0, a*]. 
Proof. We first observe that the function a H> LRL(a, /3)— RSL(a, /3) is decreasing. This follows immediately 



from Lemmas 20 and 
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We also claim that the function j3 1— > LRL(a, /3) — RSL(a, /3) is increasing for 
(3 ^ 37r/2. Together, these facts prove the lemma: consider two values a\ < ^ a*. Since f3 LR (a) is 
a decreasing function, we have /3 LR (ai) > /3l R (c* 2 ) ^ 2-7T — a* > 37r/2, and so 



LRL(a 1)( S LR (ai)) - RSL(ai,/3 LR (ai)) ^ LRL(ai, / 8 LR (a 2 )) - RSL(a x , /3 LR (a 2 )) 

^ LRL(a 2 , f3 LR (a 2 )) - RSL(a 2 , /3 LR (a 2 )). 



It remains to show that for (a,/3) € -B with /3 ^ the function /3 i-> LRL(a, /3) — RSL(a, /3) is 
increasing. Since we are in case B A , by Lemma [7] the point S lies on the arc SqS^ of Lg, while F lies on 
the arc F^F^ oi L F (see Figure 12(a) ). It follows that 7 L > j3 ^ 37r/2, and so cos7 L > 0. This implies 



that RSL(a, /3) = 1 — cos7 L < 1. On the other hand, by Lemma 
the claim follows. 
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LRL(a, j3) > 1 for p > 4f , and 

□ 



By Lemmas 22 and 24 if b(g?) > d+27r, then B(d) = B(d, a B , /3 LR (a B )), where either a B = or a B = a* 
or a B € (0, a*) is the unique value where LRL(a B , /3 LB (a B )) = RSL*(a B , /3 LR (a B )). The first case a B = 
can never occur, as a simple geometric argument similar to Lemma [4] shows that rsl(0, /3 lr (0)) ^ d + 27r. 
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that for 



The second case a B — a* holds only for d = y/2, as we have seen in the proof of Lemma 
d > y/2 we have LRL(a*, 2ir — a*)< 2a* + 2n = RSL* (a*, 2ir — a*), and lrl does not have a maximum at 
this point. For d = v2, however, lrl(«*, 2tt — a*) = 2a* + 2n = RSL* (a*, 2n — a*), so we have shown: 



Lemma 25. For < d < 2 if dub 3(d) > 2ir, then dubs(eQ = B(d,a B , f3 LR (a B )) 
unique value in [0,a*] where RSL*(<i, a B , f3 LR (a B )) = LRh(d, a B , f3 LR (a B )) . 



d, where a B is the 



It remains to argue the monotonicity of d 1— > duhg(d). 
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Lemma 26. For \/2 d < 2, dubs(d) decreases monotonically from dubs(-\/2) = 5n/2 — \/2 to 
dub B (d*) = 2tt for d* 1.5874. We have dub B (d) < 2tt /or d ^ d* . 

Proof. Assume the function is not monotone, so there is y2 ^ di < d± < 2 such that 2ir < dubs(di) < 
dubs(c?4). Since dubs(d) is continuous, it assumes its maximum D := max^g^^] dubs(d) on the closed 
interval [d^dj. The set {d G [d^dj | dubs(d) = D} is compact, and so assumes its infimum, say at d 3 . 
So we have dubB(d 3 ) > dub#(di) > 27r, and dubs(d) < dubs(d 3 ) for di ^ d < d 3 . 

By Lemma [25} we have dubs(d 3 ) = B(d 3 , a B , /3 B ) — d 3 , where < a B ^ arcsin(d 3 /2) and (3 B = Plr(c( b ). 
Let us define di = max{di, 2 sin(27r — /3 B )}. Then d\ ^ di < d 3 , and (a B ,/3 B ) G B A for d = d2- Then 
there is an RSL- or LRL-paih from (0, 0, a B ) to (d2, 0, /3 B ) of length at most dubs^) + d2- By Lemmas [T] 
and|8j there is then a path of length dubs^) + d from (0, 0, a B ) to (d, 0, f3 B ) for all d ^ d2, implying 
that £(d, a B ,/3 B ) ^ dubs(d 2 ) + d, a contradiction to the assumption that dub^(d 3 ) > dubs(d 2 ). 

The approximation to the value of d* has been computed numerically. For a given d, we first 



approximate a B numerically by binary search on the interval [0,a*] using Lemma 24 We can then 



compute /3 LR (a B ) and B(d). □ 

8 The Dubins cost function 

We now put all the pieces together. 

Lemma 27. For < d < \/2, dub(d) = dub^(d). The function decreases monotonically from dub(0) = 
7tt/3 to dub( v / 2) = 5tt/2 - a/2- 



Proof. By Lemma 18 dubyi(d) = A(d) — d for < d < \[2. By Lemma 
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this function decreases 



monotonically from dub^(0) = 7ir/3 to dubA(\/2) = 5n/2 — \/2, and so dubyt(d) > 2n = dubc(d). 



By Lemma 20, for any (a,/3) € B A there is an i?5L-path with length at most RSL*(a*,27r 



a*) = RLR(a*,27r — a*). In the proof of Lemma 14 we observed that for < d < v2 ; we have 
RLR(d, a* , 2ir - a*) < LRL(d, a* , 2tt - a*) < A(d) = dub^d) + d. It follows that dub B (d) < dub^(d), 
and so indeed dub(d) = dub^(d) for < d < \/2. □ 

Lemma 28. For \[2 ^ d < d* , dub(d) = dubs(d). The function decreases monotonically from 
dub B (a/2) = 5tt/2 - x/2 to dub(d*) = 2ir. For d* < d < 2, dub(d) = 2ir. 



Proof. By Lemma 19 we have dub^(d) < max{2-7r, dubs(d)} = max{dubc(d), dubs(d)}, and so dub(d) 



max{27r, dubs(d)}. By Lemma 26 we have dubs(d) > 2tt for y/2 ^ d < d*, and the lemma follows. □ 



Lemma 29. For d ^ 2, we have dub(d) = 2ir. 

Proof. Let d ^ 2 and (a,B) G A. We need to show that g(d,a,/3) < 2vr + d. If (a,j3) G C A , then 
this follows from Lemma|4j Otherwise we must have (a,f3) G B A , and we have S > tt/2 by Lemma[3] 
We choose di < 2 such that d\ > 2sin(7r — S) and di > d* and consider the configuration for 
distance d±. Since case A occurs only within the 5-range arcsin(di/2) ^ 6 ^ tt — arcsin(d!/2), we must 
be in either case B or case C, so there is a path of type RSR, LSR, RSL, or LRL of length at most 
di +dub(di) = di + 2-7T from (0,0, a) to (di,0,/3). By Lemmas [I] and [8] there is then a path from (0,0, a) 
to (d, 0, /?) of length at most di + 2ir + (d - di) = d + 2ir. □ 

We summarize some of our results in the following theorem: 

Theorem 1. The function dub(d) has two breakpoints at \[2 and d* « 1.5874. For d < \[2, dub(d) = 
dub A (d) < dub A (0) = f. ForV2^d<d*, dub(d) = dub B (d) < dub B (\/2) = \ - y/2. For d ^ d* , 
we have dub(d) = 2n. 
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